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Abstract 


The variational perturbation theory for wave functions, which has been shown 
to work well for bound states of the anharmonic oscillator, is applied to res¬ 
onance states of the anharmonic oscillator with negative coupling constant. 
We obtain uniformly accurate wave functions starting from the bound states. 
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Perturbation for wave functions of quantum systems is more complicated problem than 
that for energy eigenvalues since the former deals with functions while the latter only deals 
with a number. Even for energy eigenvalues, however, naive perturbations by a coupling 
constant give rise to divergent series in most cases and numerous optimization techniques 
have been studied |^. Variational perturbation method, sometimes called the 6 expansion, 
is one of the successful approaches. It has been found to reproduce energy eigenvalues of 
quantum anharmonic oscillator (AHO) and double-well potential (DWP) with high accuracy 
PJ. Characteristic of this method is that one introduces artihcial parameters which are not 
contained in an original Hamiltonian and determines those parameters after calculating a 
physical quantity perturbatively. One of the advantage of the method is such flexibility that 
one can choose different value of the parameter depending on the physical quantity one is 
interested in. One can also choose different value depending on the order of the perturbation 
one employs. In our previous work, we showed that this advantage could be suitably taken 
to apply the method to wave functions and obtained uniformly accurate wave functions of 
AHO and DWP 0. 

Recently, it has been shown that a variational perturbation approach can be also adapted 
to treat systems which accompany quantum tunneling decay processes 0. In these systems, 
energy eigenvalues have imaginary parts and naive perturbations fail in general. In this 
article, motivated by these studies, we attempt to apply our method in 0 to wave functions 
of a tunneling system. As a simple model, we take anharmonic oscillator with negative 
coupling constant (NAHO). 

We hrst recall one dimensional AHO to illustrate the variational perturbation treatment 
for wave functions. The Hamiltonian is 


One dehnes a new Hamiltonian A) = iPo(O) -|- 6Hi{Q, A) with 

Ho{n) = ^ (-^ + A) = i ((1 - + Xx^) ( 2 ) 

such that Hs=o{fl,X) = Hq{Q) and Hs=i{fl,X) = H{X). At hrst, one employs Rayleigh- 
Schrodinger (RS) perturbation with 6 as an expansion parameter. Bases of the expansion 
are the eigenvectors of Hq{Q) which is harmonic oscillator (HO) with a trial frequency D. 
The u-th order perturbative wave functions thus obtained have the following form 

n 

j 1=0 


Next step is to choose the parameter D so as to make approximate functions uniformly valid. 
This can be achieved 0 by 


D) 


dQ 


5=1 


0 . 


( 4 ) 


This condition can be solved for D, in principle, as a function of x (and A). By substituting 
the resultant 0(x) to (|^), we get optimized wave functions ijj^'^'^{x;^l{x)). 
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We can interpret the condition (^) in two different ways. First interpretation is a geomet¬ 
rical one. If we regard the obtained functions f2) as a family of curves parametrized 

by Eq.(^ is nothing but the necessary condition for to be an envelope of 

the family 0. The other interpretation is the principle of minimum sensitivity (PMS) P|; 
perturbatively calculated quantity should not depend on the artificial parameter at 5 = 1 
since the exact one does not depend on at all. 

For the ground state wave function, first order variational perturbation results in 

4^\x-,n)=7r-^/^n^/^exp 

(5) 

The condition (§) applied to (H) reads 

IQx^n^ + - A{2\x^ + 4x^ + 7)Q^ - 4(9Ax^ + 4x^)Q^ - 2(33Ax^ - 14)fl + 99A = 0. (6) 


X 


1 + 5 



A 


n-^){2nx^ - 1 ) — n-\ 4 n^x'^ + i2nx^ - 9 ) 
32 


Since this is a higher degree algebraic equation for one cannot solve it analytically in 
general. However asymptotic behavior of the physical solution at long distance can be 
estimated by Eq.(|^) and found to be 

> ^x"^ (x —>• cx)). (7) 

Therefore, the perturbative wave function improved by the condition (^) has its asymptotic 
form 




exp 




( 8 ) 


The result obtained above has correct x dependence and only the numerical factor 8 deviates 
a bit from the exact value 9. A noticeable success of the method in this case can be 
appreciated by this correct asymptotic behavior. One can never obtain such a behavior by 
naive RS perturbation since perturbative wave functions obtained by RS at any finite order 
result in finite superposition of the wave functions of HO. For example, the ground state 
wave function obtained by naive RS perturbation is of the following expression]^ 




(9) 


In the case of A > 0 the physical solutions of Eq. (j^) are to be restricted to real number 
but this is not the case for A < 0. The situation is observed by the asymptotic solutions ( 0 ); 
by Eq.(l^) with A < 0 we have 


0(x) 



( 10 ) 


^Note that this form is independent of whether A > 0 or A < 0. 
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that is, pure imaginary at infinity. Substituted them into Eq.(^ reads 

( 11 ) 

This oscillatory behavior is of course correct one for NAHO except for the numerical coeffi¬ 
cient. 

From the consideration above, we may expect to get resonance states of NAHO pertur- 
batively from the bound-state wave functions of HO if we make a novel generalization of 
our method to allow the trial frequency O to take complex values. In fact, for the ground 
state energy of NAHO, it was shown by Karrlein and Kleinert |@ that the imaginary part 
of the energy can be extracted precisely by this complex-valued optimization procedure for 
energy eigenvalues. 

In Fig.|^ we show as an example the ground state wave functions obtained in this way 
for A = —0.2. There are two wave functions which are complex conjugate each other. Thus, 
there appear two thick dashed lines in Fig.|^ which correspond to the imaginary parts of 
these two wave functions. The real part contributes high probability density in the vicinity 
of the potential valley and decreases as the potential barrier becomes high. Outside the 
barrier, both the real part and imaginary part oscillate as Fq.(^) indicates. Therefore, we 
successfully get qualitative behavior of the resonant wave functions of NAHO even at the 
first order perturbation. 

When we try to compare the perturbatively obtained results with the exact ones, we 
encounter two difficulties. Firstly, it is difficult to compute numerically an eigenstate which 
oscillates rapidly at boundaries. Secondly, the wave functions under consideration are not 
normalizable and therefore it becomes meaningless to compare numerical values of the wave 
functions at each point. To avoid these problems, we invoke the complex-coordinate method 
(CCM) 1^. This tool is often used for computing complex eigenvalues for resonance states 
numerically. The essence of CCM is as follows. The ordinary Schrodinger equations in 
configuration space 

Hipi^x) = Eil){x) (12) 

are defined on the real x axis. Then, one rotates this support of 'ip{x) in complex x plane 
as a: —>■ xe*®. One determines the argument of the rotation 6 such that the rotated wave 
functions 'ifj^xe''^) = U{9)'ip{x) tend to 0 at infinity, that is, 

U{9)HU-\e)'il;{xE^) = Ei){xE% (13) 



'ip{xe^^) —> 0 (|x| —>• cx)) 


(14) 


where U{9) denotes the rotation operator in complex x plane. Thanks to ([T^), numerical 
calculation turns out to be easier and the rotated wave functions become normalizable. 

In order that the wave functions of quartic AHO satisfy the boundary condition (0, 
the following inequality must holds [H 


(argix) -F -(argA) 
0 


< 


71 


( 15 ) 
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( 16 ) 


For NAHO, argA = vr, and the inequality reads 

--<e<o. 
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We fix the argument at 6^ = —tt/G for our calculations. The normalized ground-state wave 
function in CCM for A = —0.2 calculated numerically is shown in Fig.^]. 

Finally, we combine perturbation theory with CCM. The procedure is that we rotates 
Eq.@ in the complex x plane by the same argument as that for the exact one Eq.(|^ 
(6* = —7r/6 in our case). 

(17) 

The optimization condition Eq. thus changes in the same manner. If the physical solutions 
of Eq.(D are written formally as = f2(x), the solutions in the rotated coordinate are 
= Q{xe^^). Therefore we eventually get the optimized wave functions combined with 
CCM 


0^ (x) = n{xe^^)). (18) 

In Fig.|^, we show the Ist-order ground-state wave function for NAHO with A = —0.2 and 
6 = —ir/G, obtained by Eq.([l8Dp|. The comparison of it to the exact one in Fig.^ exhibits 
excellent agreement in the whole of a domain caluculated. To make more quantitative 
comparison, (squared) norm defined by 

ll^exact _ ^(1)||2 ^ J dx\(j)^^^^\x) - (19) 

are represented in Table | and the expectation value of the (rotated) Hamiltonian are rep¬ 
resented in Table The agreement is excellent as a first-order perturbation. 

For the purpose of comparing our variational results with naive RS ones, we also show the 
results obtained from RS perturbation. We calculate RS wave functions Eq.(|]) on complex 
coordinate 


( 20 ) 

with 6 = —ir/Q. In Fig.^, thus obtained with A = —0.2 is shown. Squared norm Eq.(|^ 
and the expectation value of the (rotated) Hamiltonian for are also shown in Table | 
and respectively. 

Some remarks are in order here. 

1. For A < —0.12 the physical solutions of Eq.(^ with real x can be obtained as a complex 
conjugate pair in the whole x domain calculated. They correspond to out-going and in¬ 
coming states respectively. For weak negative coupling A > —0.12, that is, high potential 


^After we calculate the optimized wave function Eq.(15), we make a (global) phase transformation 
so that the phase of the optimized wave function becomes the same as that of the exact wave 
function at the origin. 
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barrier case, there appear domains where the physical solutions split into two different real 
values and one cannot judge which solution should be taken. Therefore naive application 
of the optimization condition (^) do not work well in high-barrier case. This situation is 
similar to that for the same optimization scheme for eigenvalues of NAHO 0]. 

2. In the case of the optimized perturbation combined with CCM, the physical solutions of 
Eq. with complex x can be determined uniquely as a complex number since the coefficients 
of the algebraic equation are now complex. For A < —0.12, imaginary parts of those 
physical solutions are all positive in the whole domain calculated. For A > —0.12, on the 
other hand, imaginary parts of the physical solutions turn to be negative in a domain |x| < Xq 
with a certain point Xq. However, no difficulties arise in the method with CCM in contrast 
to in the case of the method with real x. 

3. From Table | we can see that the errors for our method take its maximum around 
A = —0.2. Even at this point, however, the error defined by Eq. ([T^) is about 0.0016, which 
indicates that our method is applicable to both the tunnel regime and the sliding regime of 
the coupling constant spaceQ On the contrary, in the case of the optimization for eigenvalues 
of NAHO, it was proved that there exists lower bound Aq so that optimized perturbation 
series converge to the exact ones for |A| > Aq 0. Therefore, it may be said that convergence 
property of the method for eigenfunctions with CCM is different from that for eigenvalues. 

4. We have also calculated the 2nd-order variational perturbation. The results only ex¬ 
hibit slight quantitative improvement from the Ist-order ones. This feature indicates rapid 
convergence property of the method although whether perturbative wave functions really 
converge into the exact ones or not is of course another problem. 

5. In view of both perturbation theory and physical problem, Zeeman and Stark effects are 
interesting problems^. These systems accompany quantum tunneling and thus the method 
presented in this paper will be applicable. 

In summary, we have generalized the variational perturbation method for wave functions 
of bound states to treat those of resonant states and applied it for NAHO. We can obtain 
qualitative nature of resonance states such as oscillatory asymptotic behavior and reproduc¬ 
tion of complex energy, only by a Ist-order perturbation based on the bound states of HO. 
The method with CCM can overcome the difficulty in the method with real x for A > —0.12 
and can give uniformly accurate wave functions in both the tunneling and the sliding regime. 
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^ Some authors characterize the tunneling (sliding) regime of the coupling constant space as the 
real part of the ground state energy lower (higher) than the potential barrier height and we follow 
this definition here. For our Hamiltonian Eq.(||), the border of the two regime is at A ~ —0.3387. 

^An interesting resummation approach for Stark effect is found in Ref. |^. 


6 



REFERENCES 


[1] Large Order Behavior of Perturbation Theory, Current Physics - Sources and Comments, 
vol.7, ed. J. C. Le Guillou and J. Zinn-Justin (North-Holland, Amsterdam, 1990). 

[2] See references quoted in G. A. Arteca, F. M. Fernandez and E. A. Castro, Large Order 
Perturbation Theory and Summation Methods in Quantum Mechanics, (Springer-Verlag, 
Berlin, 1990) Chap. 5, and H.Kleinert, Path Integrals in Quantum Mechanics, Statistical 
and Polymer Physics, 2nd. edition (World Scientihc, Singapore, 1995) Chap. 5. 

[3] T. Hatsuda, T .Kunihiro and T. Tanaka, Phys. Rev. Lett. 78 (1997) 3229. 

[4] H. Kleinert, Phys. Lett. B300 (1993) 261; R. Karrlein and H. Kleinert, Phys. Lett. A187 
(1994) 133. 

[5] P. M. Stevenson, Phys. Rev. D23 (1981) 2916. 

[6] C. M. Bender and T. T. Wu, Phys. Rev. 184 (1969) 1231. 

[7] See references quoted in e.g. Y. K. Ho, Phys. Rep. 99 (1983) 1, and Lecture Notes in 
Physics, vol.211, ed. S. Albeverio, L. S. Ferreira and L. Streit (Springer, Berlin, 1984) 
p.235. 

[8] R. Guida, K. Konishi and H. Suzuki, Ann. Phys. 241 (1995) 152; ibid. 241 (1996) 109, 
H. Kleinert and W. Janke, Phys. Lett. A206 (1995) 283. 

[9] K. Hiraizumi, Y. Ohshima and H. Suzuki, Phys.Lett. A216 (1996) 117. 


7 



FIGURES 


FIG. 1. The perturbative wave functions {x, for the ground state with A = —0.2. 

The thick solid line is their real part and the thick dashed lines are their imaginary parts. (See 
text for the meaning of two lines.) The thin dashed line is the potential with A = —0.2. Relation 
between types of lines and their meanings is common to all the figures in this article. 


FIG. 2. The exact wave function for the ground state with A = —0.2 

calculated by GGM. The rotation argument 9 equals —tt/G. 

FIG. 3. The variational perturbation wave function 4>Y^p{x) combined with COM for the ground 
state with A = —0.2. 


FIG. 4. The RS perturbation wave function combined with COM for the ground state 

with A = —0.2. 



TABLES 



||^exact_0W ||2 

||0exact _ 0(1) ||2 

A = -0.05 

1.29994 E-5 

1.04105 E-3 

A = -0.06 

3.34356 E-5 

1.98414 E-3 

A = -0.07 

7.73821 E-5 

3.38198 E-3 

A = -0.08 

1.59194 E-4 

5.31350 E-3 

A = -0.09 

2.89919 E-4 

7.84056 E-3 

A = -0.1 

4.75445 E-4 

1.10032 E-2 

A = -0.2 

1.56666 E-3 

7.65252 E-2 

A = -0.3 

1.18456 E-3 

1.90267 E-1 

A = -0.4 

9.16755 E-4 

3.32428 E-1 

A = -0.5 

7.57480 E-4 

4.84690 E-1 

A = -0.6 

6.56376 E-4 

6.33762 E-1 

A = -0.7 

5.87611 E-4 

7.72179 E-1 

A = -0.8 

5.38161 E-4 

8.96876 E-1 

A = -0.9 

5.01007 E-4 

1.00742 EO 

A = -1.0 

4.72111 E-4 

1.10469 EO 


TABLE I. Squared norm defined by Eq.(19). Since both 
unity, the squared norm takes a value in the interval [0,2], 


^exact 0(1) 


are normalized to 
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exact 

{H)^vp 

m'il 

A = -0.05 

0.47912 + i 0.00001 

0.47971 + i 0.00036 

0.50328 - i 0.00369 

A = -0.06 

0.47416 + i 0.00006 

0.47519 + i 0.00061 

0.51141 -i 0.00176 

A = -0.07 

0.46879 + i 0.00026 

0.47050 + i 0.00098 

0.52287 + i 0.00248 

A = -0.08 

0.46297 + i 0.00077 

0.46564 + i 0.00150 

0.53799 + i 0.00961 

A = -0.09 

0.45677 + i 0.00176 

0.46062 + i 0.00219 

0.55712 + i 0.02025 

A = -0.1 

0.45034 + i 0.00335 

0.45551 + i 0.00308 

0.58058 + i 0.03497 

A = -0.2 

0.39744 + i 0.04471 

0.39512 + i 0.03694 

1.11140 + i 0.52301 

A = -0.3 

0.37377 + i 0.09506 

0.36928 + i 0.09092 

2.28179 + i 1.88598 

A = -0.4 

0.36441 + i 0.13868 

0.35999 + i 0.13638 

4.06817 + i 4.22403 

A = -0.5 

0.36146 + i 0.17576 

0.35740 + i 0.17444 

6.30746 + i 7.37246 

A = -0.6 

0.36169 + i 0.20780 

0.35798 + i 0.20704 

8.81116 + i 11.06949 

A = -0.7 

0.36363 + i 0.23600 

0.36022 + i 0.23560 

11.43170 + i 15.07917 

A = -0.8 

0.36655 + i 0.26122 

0.36339 + i 0.26106 

14.07565 + i 19.23405 

A = -0.9 

0.37004 + i 0.28407 

0.36710 + i 0.28408 

16.69183 + i 23.42979 

A = -1.0 

0.37387 + i 0.30499 

0.37111 + i 0.30512 

19.25816 + i 27.61088 


TABLE II. Expectation value of (rotated) Hamiltonian with respect to the exact wave function 
(H)ex3-ct that with respect to the perturbative one 
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